JOURNAL CF RESEARCH of the National Bureau of Standards— D. Radio Propagation 
Vol. 67D, No. 5, September-October 1963 

Reflection of VLF Radio Waves From an Inhomogeneous 
Ionosphere. 1 Part II Perturbed Exponential Model 

James R. Wait and Lillie C. Walters 



Contribution From Central Radio Propagation Laboratory, National Bureau of Standards, Boulder, Colo. 

(Received April 10, 1963) 

The oblique reflection of radio waves from a continuously stratified ionized medium 
is considered. In this paper the medium is assumed to be isotropic. The height profile of 
the effective conductivity is a Gaussian curve superimposed on the (undisturbed) exponential 
form. The reflection coefficient is shown to be influenced by the vertical location of the 
Gaussian perturbation. In some cases the magnitude of the reflection coefficient is increased 
while, in other situations, it is decreased. In nearly all cases, insofar as phase is concerned, 
the presence of the perturbation corresponds to a lowering of the reflection height. 



1. Introduction 

In a previous communication from the present 
authors [Wait and Walters, 1963], oblique reflection 
of (VLF) radio waves from a continuously stratified 
ionized medium was considered. The profile of the 
effective conductivity was taken to be exponential 
in form. Actually, this is a fairly good representation 
of the actual D layer of the ionosphere under day- 
time conditions. Henceforth, that paper will be 
referred to simply as (I). 

It is the purpose of the present paper to consider 
profiles which are no longer exponential in form. 
Since the objective is to gain insight into the mecha- 
nism of reflection from perturbed layers, a number 
of idealizations are made. First, it is assumed, 
under quiescent conditions, that the ionospheric 
conductivity varies exponentially with height. Then 
the idealized perturbation is assumed to have a 
Gaussian form. Again, for sake of simplicity, the 
earth's magnetic field is neglected as in (I). This is 
well justified when considering effects which result 
from ionization in the lowest ionosphere. 

2. Description of the Profile 

The notation follows that used in (I) as closely as 
possible. Thus the undisturbed profile, as a func- 
tion of height z, is defined by the conductivity 
parameter 1/L(z) where 



11 r*\ 

— =I exp(^), 



(1) 



and L is a constant, is a gradient parameter and z 
is the height above the reference level z=0. Under 
the isotropic assumption, it is known that [Wait, 
1962] 

= co(H^) ; 



L= 
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in terms of the angular frequency co, collision fre- 
quency v> and plasma frequency oo . At VLF, 
^>>co, and therefore 



T CO , C0 

L^ — where co r = — ; 

C0 r V 



(3) 



to within a very good approximation. 



In general, it is seen that 



Liz) 



is proportional to 



where N(z) and v(z) are the electron density 



N(z) 

v{z) 

and collision frequency regarded as a function of 
height. The constant 0, in the exponent, is a meas- 
ure of the sharpness of the gradient. For example, 
0=1 km -1 means that the ratio of co r (2) or N(z)jv(z) 
increases by 2.71 for each km of vertical height. 
From the recent work of Barrington et al. [1962], 
Kane [1962], and Belrose [1963], it appears that /3 
for an undisturbed ionosphere may be in the range 
from 0.2 to 0.8. If the level from about 60 km to 
70 km is considered, it appears that 0=0.3 typifies 
many of these daytime ZMayer profiles. A detailed 
study of the influence of changing is to be found 
in (I). For this paper, /3 is chosen to be 0.3. 

Having specified our undisturbed profile, we now 
wish to introduce the perturbation. It is assumed 
that the collision frequency profile is unchanged 
whereas the ionization is to be increased by an 
amount AN(z) where 



AN(z)=AN exv 



n^n- 



(4) 



and AiVo, F and D are constants. Clearly, the 
maximum value of AN(z) is AN Q which is located 
at z=F. Furthermore, the thickness of this layer 
is 2D which is the vertical distance between the 
levels where AN(z) drops to AN /e. 

In order to estimate correctly the influence of this 
Gaussian shaped layer, it is necessary to assume 



519 



something about the collision frequency profile. A 
careful study of the recent literature indicates that 
an exponential variation of v(z) with height z is not 
unreasonable. The form chosen here is 



v(z)=v exp (- 



(5) 



where /3=0.3 Ion l . Therefore, the resulting con- 
ductivity perturbation has the form 



AN(z) = AN 
v{z) ' 



vo 
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Figure 1. The undisturbed and disturbed conductivity profiles 
used in this paper. 



The complete profile, under these idealized con- 
ditions, is given by 



+A exp (J z) exp [-(£Z_T)^ 



(7) 



where the right-hand side is proportional to the 
effective conductivity of the medium as a function 
of height above (or below) the reference level at 
2=0. The coefficient A defines the strength of the 
perturbation. In fact, 



A- 



A,Vo 



where N is the electron density of the undisturbed 
profile at the reference level z=0. In this paper, 
as in (I), Z(0) = 7.5/\ where X is the wavelength in 
kilometers. 

It is admitted that other ways to define a pertur- 
bation in the profile may be preferable. Here the 
electron density anomaly, for a given value of A, 
does not change with its vertical location F. Con- 
sequently, we may anticipate that the influence of 
this type of perturbation will be diminished at 
sufficiently low heights because of the increasing 
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Figures 2a and b. The reflection coefficient as a function of the vertical location, F, of the Gaussian perturbation, for various 

angles of incidence. 
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collision frequency. However, we shall see the 
problem is not quite this simple as other factors 
come into play. 

A sketch of the profiles used is given in figure 1. 
The undisturbed profile is the exponential form, 
while the disturbed profiles have the superimposed 
Gaussian "bump." The location of the "bump" 
for five typical profiles is specified by the appropriate 
value of F. 

3. Results of the Calculations 

The method used to calculate the reflection co- 
efficient R has been described in detail in (I). The 
quantities considered are the amplitude \R\ and the 
phase of R for a vertically polarized plane wave 
incident at an angle whose cosine is C. The reflec- 
tion coefficient is evaluated in the free space region 
corresponding to z->— °o . However, it is important 
to remember that the phase is referred to the level 
z=0. 

The plan of the calculations is to vary the value 
of one parameter while keeping the others constant. 
To obtain a complete understanding of the various 
phenomena, an enormous number of calculations is 
needed. In order to keep the problem within reason- 



able bounds and to reduce the expense of the compu- 
tation, only a limited number of cases was considered. 
These results are shown in graphical form in figures 
2 to 6. In all cases 13=0.3 km -1 . 

In figure 2a the amplitude of the reflection co- 
efficient is plotted as a function of F for \=15 km 
(/=20 kc/s), A = 2, D=2 km, and (7 values varying 
from 0.05 to 0.4. Small values of (7 here correspond 
to angles near grazing. For long distance propaga- 
tion of VLF radio waves, values of G near 0.1 are 
most important. For this case, it is interesting to 
note, when F is near or above zero, that \R\ takes 
the same value as for the undisturbed profile. As 
the "bump" or perturbed layer is lowered, the reflec- 
tion coefficient first increases then decreases. Even- 
tually, as the "bump" is brought down to very low 
heights, \R\ returns to its undisturbed value. The 
other curves for highly oblique incidence have a 
similar behavior. Thus, the "bump" may either im- 
prove or degrade the reflection. Presumably, at the 
lower heights the Gaussian layer is acting as an 
absorber whereas, at greater heights, it enhances the 
reflection. At the steeper angles of incidence, the 
situation becomes more complicated. It is probable 
that this results from interference between multiple 
reflected rays between the upper side of the "bump" 
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Figures 3a and b. The reflection coefficient, as a function of F, for various widths of the Gaussian pertubation when C = 0.2 {i.e., 

angle of incidence is 78°). 
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Figures 4a and b. The reflection coefficient, as a function of F, for various widths of the Gaussian perturbation when C = 0.1 (i.e., 

angle of incidence is 84°). 
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Figures 5a and b. The reflection coefficient, as a function of F, for various wavelengths (from 10 to SO km). 



and the exponentially varying layer. Such an inter- 
ference phenomenon becomes more pronounced at 
steeper incidence because the vertical component of 
the wavelength is becoming comparable with typical 
values of F. 

The phase of R is shown in figure 2b for the same 
conditions as in figure 2a. Again, it is apparent that, 
when the Gaussian "bump" is such that F is near 
or above, the phase of R attains its undisturbed value. 
When the angle of incidence is highly oblique the 
phase undergoes an increase (i.e., decrease of lag) 



as the "bump" comes down to lower heights. Suffi- 
ciently far below the reference level, the phase of R 
returns to its undisturbed value. It is well to note 
that as C becomes small (i.e., approaching grazing 
incidence), the phase of R is approaching — 180°. 

For highly oblique incidence the influence of the 
"bump" is to lower the effective height of reflection 
for the whole range of F. However, a very interest- 
ing phenomenon occurs at steeper incidence. As can 
be seen in figure 2b, when (7=0.2 the phase undergoes 
a rather rapid change as F varies from about —22 
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Figures 6a and 1). The reflection coefficient, as a function of F, for various amplitudes A of the perturbation. 



km to —27 km. As C is increased further there is 
an apparent discontinuity when the phase changes by 
360°. Such a change of 2tt radians is quite permis- 
sible since the ordinate is arbitrary to within any 
integral number of 2w radians. Thus, the phase 
curves for (7=0.3 and 0.4 could have been drawn 
in the range below —160°. 

The curves in figure 2b, even if they show nothing 
else, demonstrate that phase shifts in reflection 
phenomena may have some unusual cycle am- 
biguities. 

The influence of the width of the Gaussian pertur- 
bation or "bump" is shown in figure 3a for the 
amplitude 17*1 and in figure 3b for the phase of E. 
Here A = 2, X=15 km, and (7=0.2. The amplitude 
curves show that, when I) is increased, the overall 
influence of the layer becomes somewhat greater. 
There is some tendency for the thinner layers (i.e., 
smaller D) to be more effective at greater heights. 
The corresponding phase curves show that the 
thicker layers always produce a larger phase change. 
Furthermore, as D exceeds 2 km, a point is reached 
where the "360° jump" takes place. 

The curves in figures 4a and b are for the same 
conditions as figures 3a and b except that now C= 
0.1, corresponding to nearer grazing incidence. 
The amplitude curves have a very similar shape. 
The phase curves are also similar except that the 
"360° jump" is no longer present. 

The wavelength dependence of the reflection 
coefficient is shown in figures 5a and b. For these 
(7=0.1, A=2, and Z>=2. The wavelengths chosen 
(10, 15, 20, 25, 30km) correspond to frequencies of 
30, 20, 15, 12, and 10 kc/s. Qualitatively, the 
curves have a very similar shape. There is some 
tendency for the shorter wavelengths to be accom- 
panied by more pronounced changes. In all cases 
the "bump" acts as an absorber at low heights 
while it enhances the reflection at greater heights. 



Finally, in figures 6a and b, the influence of A, 
the relative magnitude of the anomalous electron 
density, is shown. As expected, the individual 
curves are similar in shape with the larger values of 
A corresponding to an increased change over the 
undisturbed values. It is important to note thai 
the phase anomaly is almost directly proportional 
to A. 

4. Final Remarks 

The results given here constitute a small portion 
of extensive computations dealing with reflection 
of waves from inhomogeneous media. In subse- 
quent parts to this series other types of profiles will 
be considered. Also, the applications to the mode 
theory of VLF propagation are to be described in 
some detail. 
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